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Abstract 
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1 Introduction 



Glassy dynamics is often described as motion in a complex energy landscape M, with 
many valleys separated by barriers. Valleys have been given various definitions [0, [| f|. |5J . 
which are not equivalent from a dynamical point of view pj. One common feature of all 
these approaches is the exponential growth with the system size of the number of valleys 
at fixed energy density E: 

N(N;E)~exp(NS ap (E)), (1.1) 

where S ap (E) is the 'a priori' configurational entropy, or complexity. 

As valleys are so numerous, an important issue concerns the possibility of giving an a 
priori statistical description of them. The key question to address is then the following. 
Assuming the initial configuration be random, does the system sample all the possible 
valleys with equal statistical weights, i.e., with a uniform or flat measure, or, to the 
contrary, is the size of the basin of attraction of each valley relevant? The assumption 
that an a priori statistical description with a flat measure could hold was first formulated 
by Edwards |7| in the context of granular materials, and is therefore referred to as the 
Edwards hypothesis. 

One-dimensional spin systems at zero temperature provide an adequate setting for 
the investigation of the dynamical role of valleys, and especially for testing the Edwards 
hypothesis. Indeed, the blocked configurations reached by zero-temperature dynamics are 
truly metastable states with infinite lifetimes. In an earlier work jSj we have addressed 
this question for Ising chains evolving under irreversible zero-temperature dynamics, 
where each spin may flip at most once before a blocked state is reached. These dynamical 
models can be exactly mapped onto random sequential adsorption problems || , for which 
analytical tools are available. 

In the present work we pursue the investigation of metastable states in Ising chains, 
by considering the richer situation of zero-temperature Kawasaki (spin-exchange) dy- 
rP[ . We allow the rate Wq of diffusion processes (constant-energy moves) to 



namics 



vary |TT|. If the kinetic constraint Wo = is imposed [|r2j| , only irreversible processes 
are present M. For the unconstrained dynamics (Wo > 0), where diffusion processes 
are allowed, the dynamics is partly reversible. Each spin only flips a finite number of 
times, before the system globally reaches a blocked state after a finite time. The statistics 
of this blocking time is studied in Section |2|. Our main goal is then again to test the 
Edwards hypothesis, i.e., to evaluate various observables in the blocked configurations 
reached by the dynamics, starting from a random non-magnetized initial configuration. 
We shall systematically compare the data of extensive numerical simulations, presented 
in Section |], with the predictions of the a priori approach, derived in Section ^. Finally, 
some aspects of persistence, in particular the distribution of the number of flips of a given 
spin, are discussed in Section B[ 
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To be more specific, we consider a ferromagnetic chain of Ising spins a n = ±1, whose 
Hamiltonian reads 

H = -^2a n a n+1 . (1.2) 



-+), 



In Kawasaki dynamics [10|, only pairs of opposite spins may be flipped (H — * 
so that the magnetization is locally conserved. 

For simplicity, we limit ourselves to Monte-Carlo dynamics with random sequential 
updating, and we assume that the flipping rate only depends on the energy difference STC 
involved in the proposed move. Any given pair of opposite spins thus has a probability 
Wsh °f flipping per unit time, with 

5H = 2(a n _i(T n + <7 n+1 a n+2 ) e {-4, 0, 4}. 

The requirement of detailed balance with respect to the Hamiltonian ( |1.2| ) at temperature 
T = 1/(3 yields a single condition, 



W 4 

wl. 
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upon the rates W4, W , 1K_ 4 . 

We furthermore restrict ourselves to zero-temperature dynamics, hence W4 = 0. We 
choose time units such that W-4 = 1, keeping Wq as a free parameter. The allowed 
moves and the corresponding rates are listed in Table |l|. The zero-temperature limits of 
the Metropolis and heat-bath rules correspond respectively to W = 1 and Wq = 1/2. 
Hereafter we focus our attention on the range < Wq < 1. 



5H 


type 


name 


rate 


moves 


-4 


irreversible 


condensation 


1 


- + -+->-- + + 

+ - + ► + + -- 





reversible 


diffusion 


W 


+ + -+ <-> + - + + 
- + -_<-> + - 



Table 1: Allowed moves in zero-temperature Kawasaki dynamics. 



2 Statistics of blocking time 

Let us first recall that for the kinetically constrained model (Wq = 0) a finite system con- 
sisting of iV spins reaches a blocked configuration after a finite blocking time T/v ~ In TV, 
which is the jamming time of the equivalent problem of random sequential adsorption 
of hollow trimers |§ . The blocked configurations are characterized by the property that 

the spin patterns H 1 — and — I h are absent. Equivalently, there are at most two 

consecutive unsatisfied bonds. 



3 



In the present case (Wo > 0) the system still gets trapped in a blocked configuration. 
However the diffusive motion of free + spins in domains of — spins, and vice-versa, is 
allowed. Each free spin will eventually be annihilated, by meeting either another free spin 
or one of the boundaries of the domain. Blocked configurations of Kawasaki dynamics 

are therefore characterized by the property that the patterns H h and — I — are absent. 

Equivalently, isolated spins are absent, or unsatisfied bonds are isolated. 

In order to understand the statistics of the blocking time T/v, we consider first the 
regime Wo <C 1, where the time scales of condensation and diffusion are well sepa- 
rated [fLlJ| . The fast part of the dynamics, which occurs with unit rate, is identical to 
the constrained dynamics considered in Refs. || 0. For intermediate times of order 
1 <C t <C 1/Wq, the system is therefore approximately left in one of the final configura- 
tions of the constrained dynamics. The slow, diffusive part of the dynamics then takes 
place at rate W <C 1. The late stages of the dynamics are governed by large domains, 
on which a single free spin diffuses. We assume that such large domains, of size L ^> 1, 
occur with an exponentially small probability 

fdu(L) ~ exp(-L/£ dif ), 

where £dif is the relevant characteristic length. The spin diffusion constant reads D = Wq 
in our units. For a random initial point, the survival probability is known [13| to decay 
as S(t;L) « (8L/vt 2 ) exp(— 7i 2 Wot/L 2 ). The mean density of free spins at time t can 
therefore be estimated as 

S(t) « *£f m {L)S{t;L) ~ /"exp (-^- - LdL. (2.1) 

Evaluating the integral by the method of steepest descent, we thus obtain stretched 
exponential decay for the density of free spins: 

S(t) . t V» exp (-SM1) , (2.2, 



with 



A 



This behavior, already emphasized in Ref. |TT[, is a general characteristic feature of 



diffusion processes in the presence of random traps jn|. The exponent 1/3 is related to 
the one-dimensional geometry; it would read dj [d + 2) for trapping problems in higher 
spatial dimension d. The above result should hold for the late stages (Wot ^> 1) of 
Kawasaki dynamics, for any finite value of Wo- The characteristic length £dif is expected 
to depend smoothly on Wq. 

For a finite system of N spins, the last free spin will be annihilated at a time T/v such 
that NS(Tn) ~ 1, hence (WqTn) 1 ^ ~ AaiflniV. More precisely, as the histories of spins 
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diffusing on different domains are statistically independent, it can be argued along the 
lines of H that Tjv is given according to extreme-value statistics |L4 



as 



(PTOv) 1/3 « Auf (ln(iV(lniV) 3 / 2 ) + & + (2.4) 

In this expression, the factor (lniV) 3//2 takes account of the t 1 ! 2 prefactor in the survival 
probability (2.2), while the effective constant b encompasses all sub leading effects that 
have been neglected, and the random variable is distributed according to the Gumbel 
law 

/(X)=exp(-X-e- x ). (2.5) 

These predictions have been checked against extensive numerical simulations, per- 
formed according to the zero-temperature Kawasaki dynamics summarized in Table [I], 
with random sequential updates and periodic boundary conditions, starting from a ran- 
dom initial configuration with zero magnetization. The dynamics is run until the system 
gets trapped in a blocked configuration. The blocking time T/v is recorded for each sam- 
ple. This measurement has been performed for many samples (10 8 spins in total for each 
value of Wq and of N). The mean {(WqTn) 1 ' 3 ) is found to follow an almost perfect linear 
law when plotted against ln(iV(lniV) 3//2 ), at least for N ranging from 50 to 3200, whereas 
the same data plotted against IniV are bent in a significant way. Using fl2.3|) , the slope 
of the latter plot yields the value of the characteristic length £ dif , which is plotted in Fig- 
ure [l], against Wo- The error on this estimate, containing a systematic and a statistical 
part, is roughly comparable to the symbol size. The length £dif exhibits a rather weak 
dependence on Wo, decreasing from £dif = 1-77 in the Wo — > limit to £dif = 1-61 for 
W = l. 

Another confirmation of the above picture is provided by Figure |2], showing a his- 
togram plot of the variable defined by (|2.4Tj , for 125,000 samples of size N = 800, 
with Wo = 1. The parameter = 0.339, i.e., £ d if — 1-61, is taken from the data 
of Figure [I], while the constant b is chosen by fitting the average (Xn). A convincing 
agreement is found with the limit law (|2.5|). 



3 Blocked states: a priori ensemble 

As already stated, the blocked states of zero-temperature Kawasaki dynamics are the 
configurations where isolated spins are absent. Equivalently, unsatisfied bonds are iso- 
lated. 

This section is devoted to the statistical description of the blocked configurations 
taken with equal a priori weights (Edwards approach). We shall distinguish the full 
ensemble of all the blocked configurations, irrespective of their energy, the restricted 
ensemble of blocked configurations with prescribed energy density E, and the canonical 
ensemble of blocked configurations, obtained by fixing the parameter (3 conjugate to the 
energy density E (see below). 
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1.80 




Figure 1: Plot of the characteristic length £<jif> extracted from the size dependence of 
((WoT/v) 1//3 ), against the rate Wo. Symbols: numerical data. Errors are comparable to the 
symbol size. Line: third-degree polynomial fit. 



3.1 Statistics of energy and configurational entropy 



For a finite chain of N spins, we first investigate the number J\f(N; E) of blocked con- 
figurations with prescribed energy density E. This number can be evaluated by an 
elementary combinatorial reasoning [KJ. For the sake of generality, we prefer to resort to 
the transfer-matrix formalism. We introduce the partition function Z^(f3), denned as a 
sum over all the blocked configurations C = {cr n } of the Boltzmann weight associated 
with the Hamiltonian (L2). We have 



Z N {(3) 



c 



-/3H(C) 



JM(N; E)e~ mE dE. 



(3.1) 



In this framework, the parameter (3 can be positive or negative, and it is not related to 
physical temperature. The transfer matrix is then a very useful tool. Indeed the partition 
functions Zfi, labeled by the prescribed value a = (ctjv-i, ^n) of the last two spins, obey 
the recursion relation ++ ++ 



7-+ 



J N 



\z- N 
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where the 4x4 transfer matrix T reads 



r 





V 















\ 


-0 



e' 3 



The characteristic polynomial of T factors as det(Al— T) = (A 2 — e^A — 1)(A 2 — er A+l). 
The first (resp. the second) factor yields eigenvalues Ai, A2 (resp. A3, A4), with 



Al,2 



,0 



± x/e 2 / 3 + 4 



A 



3,4 



e 2,s _ 4 



The left and right eigenvectors (L a \ and are even (resp. odd) for a = 1, 2 (resp. a = 
3,4) under the spin symmetry + «-> — . We assume that they are normalized so that 
(L a \Rb) = 6 a b- Their explicit expressions will not be needed hereafter. 

For large N, we have Z N (j3) ~ A^, as Ai is the largest eigenvalue. Using ( |3.1| ), we ob- 
tain an exponential law of the form ( |1 . 1| ) for J\f(N; E), where the a priori configurational 
entropy S ap (E) is related to lnAi(/3) by a Legendre transform: 



S ap (E)- In \ 1 (P)=f3E, E 



dlnAi 

d/3 ' 



dS : 



ap 



dE 



Explicitly, we have 

E = --4=, J = -^2*=, (3.2) 
and the a priori entropy reads 

S ap (£) = E \n(-2E) + ln(l - E) - ln(l + E). 

This entropy is non-zero for — 1 < E < 0. It takes its maximal value 

S* = ln$ = 0.481212, 
where $ = (1 + VE)/2 is the golden mean, for 

E* = — \= = -0.447214, (3.3) 
V5 



corresponding to j3 = 0. Equation ( |3.3|) therefore yields the typical a priori energy density 
of a blocked configuration. 

The result ( |1 . 1|) can be recast as follows. Consider the full ensemble of all the blocked 
configurations, irrespective of their energy. The probability of observing, in that ensem- 
ble, a blocked configuration with energy density E reads 

P ap (£) ~ exp(-iV£ ap (£)), £ ap (£) = S* - S ap (E). (3.4) 

The function H ap (E) vanishes quadratically as 

5v/5 



Z ap (E)*c(E-E*y 
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The bulk of the a priori distribution of E is therefore a narrow Gaussian around E*, 
whose rescaled variance asymptotically reads 

iV Var E = N((E 2 ) - E* 2 ) — > — = — = 0.357771. 

n^oo 2c 25 

The above results also allow to determine the higher cumulants of the energy in the 
canonical ensemble (at fixed parameter f3). One has indeed 

1 , / sne\ 1 i Zn{P ~ s) Xl{(3 - s) 

By expanding this result as a power series in s, we obtain explicit expressions for the 
cumulants ((E k )) function of (3. 

As far as mean quantities are concerned, the microcanonical or restricted ensem- 
ble (fixed energy density E) and the canonical one (fixed conjugate parameter (3) are 
equivalent. It is therefore justified to recast canonical results in terms of the prescribed 



8 



value E of the mean energy, using ( |3.2|) . Generalizing this equivalence prescription to 
the cumulants of the energy, we obtain 



N Var E = N((E 2 )) — ► -E(l - E 2 ), 

N 2 ((E 3 )) — > -E(1-E 2 )(3E 2 -1), (3.5) 

iV— *oo 

N 3 ((E 4 )) — ► -E(l -E 2 )(15E 4 - 12E 2 + 1). 

iV— >oo 

The results ( |3.5|) can be given the following interpretation. Fixing (3 amounts to fixing the 
extensive part NE of the energy, while the equivalence prescription is a natural Ansatz 
to describe the fluctuations of its non-extensive part. 



3.2 Spin correlation function 



The spin correlation function C n = (<Jo<J n ) in the canonical a priori ensemble can also be 
evaluated by the transfer-matrix method. In the bulk of an infinitely long chain, and for 
n > 0, we have 



£<£i|«S|iO(£a|«S|i2i> 



A. 



(3.6) 



where 



S 



(I 



Vo 











1 





o \ 






is the spin operator. 

Because of symmetry, only the eigenvectors with a = 3, 4 contribute to the above sum. 
The values Co = 1 and C\ = —E allow to determine the products of matrix elements 
entering ( |3.6| ), without knowing the eigenvectors explicitly. The following alternative 
reasoning can also be used. Being a linear combination of {X^/Xi) 11 and (A4/Ai) n , C n can 
be shown to obey the three-term recursion relation 



1 - E)C n+ 2 + 2EC n+1 + (1 + E)C n 



0. 



The initial values Cq = 1 and C\ = —E are therefore again sufficient to determine the 
correlation function for all values of the distance n. 

For — l/v2 < E < 0, which contains the range of final energies reached by the 
dynamics, A34 = exp(±iQ) are complex numbers with unit modulus, with 



tanQ 



VI - 2E 2 



-E 



(0 < Q < tt/2). 



We are thus led to the expression 

'1 + E\ n/2 



Cm 



E 



cos nQ + 



E 2 
VI - 2E 2 



sin nQ 



(3.7) 



(3.8) 
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for n > 0. In the a priori ensemble, the spin correlation function therefore exhibits an 
exponential decay, of the form exp(— n/£ spin ), modulated by oscillations at wavevector Q. 
Both the correlation length 

£spin = i p (3-9) 

In— - 

l + E 

and the wavevector Q, given by ( |3.7| ), depend continuously on energy 

We finally quote for further reference the value of the reduced susceptibility 

oo 

X= E °n- (3-10) 

n=— oo 

After some algebra, we obtain the simple result 

^ —E(l — E) 

3.3 Distribution of domain sizes 

Another characteristic feature of blocked states is the distribution of domain sizes f(£), 
defined as the probability that a given domain consists exactly of £ consecutive parallel 
spins. Since isolated spins are absent, domains have at least size two {£ > 2). 

In the a priori ensemble, the distribution f(£) can again be evaluated by the transfer- 
matrix method. Indeed p(£), the density (per unit length) of domains consisting of exactly 
£ spins, admits an expression similar to the middle side of (|3.6|), namely 



where 



l + E 


(r 


-2E y 


1 


-E 




-EJ 




(0 





0\ 










1 







1 







Vo 





0/ 



v= (l+-)(-+l) + (l -+)(+-!) 



is the appropriate domain boundary operator. The explicit expression of the matrix 
element is not needed, as it can be fixed by normalization. The probability distribution 
f(£) of domain sizes indeed reads f(£) = p(£)/p, where p is the total density of domains 
(or equivalently, of domain walls). 

We thus obtain the geometric (i.e., discrete exponential) probability distribution 

f{i) = i±| ff!£f 2 (3,2) 



1 - E VI- E 

for £ > 2, and consistently 



00 1 i i ip 

5>0 = 7a = -V = (^l^i^i)- (3-13) 



1=2 
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The characteristic length of ( p. 12 ), 

£dom = 1 — E ' (3-14) 

In 

-2E 

is in general different from the spin correlation length £ spin of ( |3.9|) , except for the typical 
value of energy (fUSD, where £ spin = £ dom = 1/S*. 

The mean of the domain size distribution ( |3.12| ) agrees with ( |3.13| ), while its variance 
reads 

Consider now a large sample of N spins, in the canonical a priori ensemble (fixed 
parameter 0). The number M of domains in the sample is such that N = l\ + ■ ■ • + £m is 
the sum of M independent variables distributed according to f(£), neglecting boundary 
effects. The expected number of domains and its variance are given by 



(M) 1 1 + E VarM Var£ —E(l — E 2 



N n^oo (£) 2 ' N ^v-oo (£)3 4 

One has therefore 



(3.15) 



„ 4VarM 4Var£ 
NV&yE = > — — - 3.16 

N N^oo (£)3 

so that the second result of (|3.15|) agrees with expression (|3.5|) for NVarE. 



4 Blocked states: dynamics 

In this section we compare the predictions of the a priori approach, derived in Section |3], 
to the results of numerical simulations concerning the blocked configurations reached 
by the dynamics. We have used the rules of the zero-temperature Kawasaki dynamics 
summarized in Table [I], with random sequential updates, starting from a random non- 
magnetized initial configuration. 



4.1 Mean energy 

The first and the simplest quantity to be measured is the mean energy E of the blocked 
configurations reached by the dynamics. 

Let us start with a reminder of the kinetically constrained model (Wq = 0). For a 
random initial configuration, the two kinds of allowed defects, namely isolated unsatis- 
fied bonds (domain walls) and domains of two unsatisfied bond (isolated spins) occur 
respectively with the following densities in blocked configurations || [12| : 

1 / , , r 3 / 2 9 \ e" 5 / 4 

qi = -h- e ~ 5/4 - e" 9/4 J e v dy 1 = 0.219704, q 2 = = 0.071626, 
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so that the mean energy of blocked configurations reads 

E = -l + 2 qi + 4q 2 = -e" 9/4 J e y dy = -0.274087. (4.1) 

In the present case (Wo 7^ 0), the mean energy E of blocked states is expected to be 
below this number. Indeed the diffusive moves can only help relaxing more efficiently the 
energy excess of the disordered initial state. 

Figure |3| shows a plot of the mean final energy E against the diffusive rate Wq. Each 
data point corresponds to 10 8 spins in total. We have checked that no appreciable size 
dependence is to be observed. The final energy is found to be well below (|4.1| ), and below 
the typical a priori value (|3.3| ). It only exhibits a very weak dependence on the rate W , 
increasing from the extrapolated minimum value 

E(0) = -0.5279 (4.2) 

in the Wq — > limit to the maximum value 

E(l) = -0.51633 (4.3) 

for Wq = 1. 

The regime Wq <C 1 again deserves some more attention. The mean energy has a 
discontinuity at Wq = 0. Indeed the Wq — > limit of the energy, E(0) given by ( |4.2|) , 
is different from that of the constrained dynamics, Eq given by ( [4.1| ). This discontinuity 
can be analyzed as follows. In the diffusive part of the dynamics, free spins annihilate 
by meeting either each other or domain walls. The collision of a single free spin with 
a domain wall relaxes 4 units of energy, while the coalescence between n > 2 spins 
within a domain relaxes 4(n — 1) units of energy, i.e., 4(n — l)/n per spin. If all the 
free spins were annihilated in meeting domain walls, the final energy would assume the 
value E = E — 4p 2 - The efficiency of the diffusive relaxation mechanism can therefore 
be characterized by the ratio 

Eq - E(o) n n 

T) = = 1 — > — , 

4 P2 n > 2 n 

where Il n is the probability that a given free spin gets annihilated in a coalescence 
of n > 2 spins. The extrapolated value ( |4.2j ) of the final energy yields a rather high 
efficiency: r\ = 0.886. 

In view of the weak dependence of the final energy on the rate Wq (see Figure |3|), 
hereafter we restrict the numerical analysis to the case Wq = 1 of zero-temperature 
Metropolis dynamics. 

4.2 Higher cumulants of energy 

We now turn to the statistics of the energy of the blocked configurations, besides its 
mean value studied above. In analogy with the a priori estimate (|3.4|), the final energy 
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-0.515 




-0.530 




Figure 3: Plot of the mean final energy E against the rate Wo- Symbols: numerical data. 
Statistical errors are smaller than the symbol size. Line: fifth-degree polynomial fit, yielding 
the extrapolated value Q4.2Q in the W$ — > limit. 



is expected to obey a large-deviation formula of the type 

P(E) ~ exp(-iV£(£)), 

with Tj(E) being the dynamical entropy. This formula implies that the cumulants of the 
total energy scale as (((NE) k )) ~ N, just as in usual equilibrium situations. 

Instead of measuring the whole function S(-E), which would require a quite extensive 
numerical effort, we have measured the first four cumulants of the energy of the blocked 
configurations. Table ^| gives the measured values of the scaled energy cumulants for 
Wq = 1 and a random initial configuration. For comparison we also list the predictions 
of the full a priori ensemble fl3.5|) for E — E* of (|3.3| ), and of the restricted a priori 
ensemble, obtained by inserting into (|3.5| ) the observed value (|4.3| ) of the mean energy. 
The simulations have been performed on samples of various sizes ranging from N = 50 
to 200, having 10 10 spins in total. No systematic size dependence is observed. Statistical 
errors can be estimated to be of the order of one unit of the least significant digit. Both a 
priori schemes perform very unequally in predicting the energy cumulants. The variance 
is rather accurately predicted by both schemes, which perform equally poorly for the 
third cumulant, while the restricted scheme is definitely better for the fourth cumulant. 
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nnmprippl rp^nlt 
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full p\ nrinri 

1 1 ill Oj |J1 IV.y 1_ 1 


rpcf ripf pn a nrinri 

L Cull 1\j \}\Z-KX d [J J. IV.y 1_ 1 


((E)) = E 


-0.51633 


-0.44721 


-0.51633 


N((E 2 )) = NVarE 


0.3446 


0.35777 


0.37868 


N 2 ((E 3 )) 


0.031 


-0.14311 


-0.07582 


N 3 ((E*)) 


-0.46 


-0.28622 


-0.42906 



Table 2: First four scaled cumulants of the energy of the blocked configurations. Comparison 
between numerical results and predictions of the full and restricted a priori ensembles (see text). 



4.3 Spin correlation function 

We now turn to the spin correlation function C n = (cr a n ) in the blocked configurations. 
This quantity has two remarkable properties. First, since isolated spins are not allowed 
in blocked configurations, we have C\ = 1 — 2p = —E and Ci = 1 — 4p, where p is the 
density of domain walls ( |3.13| ), hence 

C 2 = 2d - 1. (4.4) 

Second, the total magnetization M = J2 n a n is exactly conserved in Kawasaki dynamics. 
Furthermore, for a homogeneous non-magnetized state, we have (M 2 ) = Nx, with the 
definition ( |3.10| ). The quantity x therefore has the same value in the initial and final 
configurations, i.e., for a random initial state, 

oo 

X= £ C n = l. (4.5) 

n=— oo 

Figure |] shows a comparison between numerical data (circles and full line), corre- 
sponding to 3 • 10 9 spins in total, and the prediction ( |3.8| ) of the restricted a priori 
ensemble at the observed mean energy ( |4.3| ) (triangles and dashed line). The a priori 
prediction and numerical data coincide both for n = 1 (by construction, as E has been 
imposed) and for n = 2 (as a consequence of Q4.4D ). For larger distances, the spin cor- 



relations oscillate in sign and fall off exponentially, in qualitative agreement with the a 
priori prediction. Both the period of oscillations and the decay length are observed to be 
slightly larger than those of the a priori ensemble. These observations remain however 
at a qualitative level. 

We also notice that the a priori ensemble fails to reproduce the identity ( |4.5| ). Indeed x 
of ( |3.11| ) equals unity for E = 1 — \/2 = —0.414214, a value which neither agrees with 



the most probable a priori energy E* of ( |3.3| ), nor with the observed mean energy 
Conversely, the prediction for the correlation function of the restricted a priori ensemble 
at energy (|4.3|) , plotted in Figure f|, has x = 1-6192, which is significantly different from 
unity. 
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Figure 4: Logarithmic plot of the absolute value of the spin correlation function |C n |, against 
distance n. Full (open) symbols show positive (negative) correlations. Circles and full line: 
numerical data. Statistical errors are smaller than the symbol size. Triangles and dashed line: 
prediction (3.8) of restricted a priori ensemble at energy (pO). 



4.4 Distribution of domain sizes 

We have determined the distribution of domain sizes f(£) by means of extensive numerical 
simulations, with a random initial state. Figure [| shows our data corresponding to 
3 ■ 10 9 spins in total. The observed distribution is not exponential, at variance with 
the prediction (532) of the a priori ensemble. The full line suggests an exponential 



asymptotic fall-off of the distribution, with a characteristic length £dom = 1.75. This 
length is close to the characteristic length of domains with a single diffusive spin, £dif = 
1.61, plotted in Figure |l[ As a matter of fact, the equality £d om = £dif is expected to 
hold in the Wq — > limit. The prediction ( ft.l2|) of the restricted a priori ensemble at 
the observed mean energy (|4.3| ) is shown as a straight dashed line, whose inverse slope 



is the prediction £ dom = 2.603 ( |3.14|) of the a priori approach. 

Another noticeable difference with the prediction of the a priori approach is that 
the sizes of successive domains in the final states of the dynamics are not independent. 
Indeed, if they were so, the relation ( |3.16| ) would hold, while the data yield 4 Var£/ (£) 3 = 
0.2790, a number significantly below the observed variance NVaxE = 0.3446, listed in 
Table ^. Hence domain sizes exhibit a weak but definite positive correlation. Domain 
sizes which are neither statistically independent nor exponentially distributed have also 
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Figure 5: Logarithmic plot of the probability distribution f(£) against domain size £. Circles: 
numerical data. Statistical errors are much smaller than the symbol size. Dashed straight line: 
prediction ( 3.1 2| ) of a priori ensemble at energy ( |4.3| ). Full straight line: guide to the eye with 
inverse slope £dom = 1-75. 



been observed recently |H| in spin chains undergoing tapping dynamics. 



5 Features of persistence 

We now turn to features related to the whole history of a single given spin, say <7q, the 
spin situated at the origin. This kind of problems belongs to the realm of persistence. In 
the present context, the central quantity is u, the total number of times <Jq flips during 
the history of the sample, before a blocked configuration is reached. The number v of 
spin flips is finite with unit probability in the limit of a large system. Zero-temperature 
Kawasaki dynamics is therefore of type JF in the classification of |Uj. Furthermore v is 
random, as it depends both on the initial configuration and on the history of the system. 
We are therefore interested in the distribution p u of the number of flips, which is expected 
to have a well-behaved limit when o"o is deep inside a large enough sample. 

Before presenting numerical data, we first predict the main salient features of the 
distribution p v , following the lines of Section |2|. Consider again a large domain of L 3> 1 
spins, containing the origin, on which a single free spin diffuses. The spin <jq flips twice 
each time the free spin traverses the origin. We are thus led to the following effective 
problem. 
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Consider a random walker in the interval —L\ < n < L 2} with absorbing boundaries 
at n = —L% and n = L 2 . The walker starts from the origin (n = 0). The probability that 
the walker returns to the origin before being absorbed by either boundary reads 

P I « l ( Ll ,L 2 ) = l-I(± + ±). (5.1) 



This formula relies on a well-known result in the gambler's ruin problem ]17|. Indeed 



suppose that the walker's first jump is to the right (resp. to the left), and consider 
< n < L 2 (resp. < —n < L\) as the gambler's wealth. Then the ruin probability 
reads P = 1 — 1/L 2 (resp. P = 1 — 1/L\). The expression ( |5.1| ) is the arithmetical mean 
of both ruin probabilities. 

The probability that a flips an even number v = 2k ^> 1 of times therefore approx- 
imately reads p 2 k(L 1 , L 2 ) ~ (P rct (L 1 , L 2 )) k {l — P ret (L l5 L 2 )). Hence the distribution of 
the number of flips can be estimated, in analogy with ( |2.1| ), as 

P2k « E ( L i + L 2 - 1) MLi + L 2 -1) P xet {L x , L 2 )\l - P rct (L u L 2 )) 
The integral entering this expression closely resembles that entering (|2.1| ). For k ^> 1 it 



is legitimate to use the steepest-descent method. The saddle point lies at L c = y/c^ di f/2. 
We thus obtain the stretched exponential law 



p 2fc ~ v/fcexp^-2^1?) (5.2) 



for the distribution of the number of spin flips, provided v = 2k is a large even number. 

The occurrence of odd numbers v = 2k + 1 of spin flips can also be explained in the 
above framework, if the spin do is either situated at an endpoint of a domain, or involved 
in a coalescence event between two free spins. Both effects are expected to scale as the 
inverse of the domain size L. This leads us to predict that odd values of v are suppressed 
by a factor of order 1/L C , i.e., 

^±i^^. (5.3) 

P2k VK^dif 

Figure ^| shows a logarithmic plot of the distribution p u of the number of flips. The 
simulations again concern samples of various sizes having 2 ■ 10 10 spins in total, with 
Wo = 1. Even numbers of spin flips v = 2k (full symbols) are clearly more frequent 
than odd numbers v = 2k + 1 (open symbols), especially for large values of k. From a 
quantitative viewpoint, the full lines on Figure ^ show a common fit of the numerical 
data for v > 10 according to the asymptotic predictions (|5.2|) , (|5.3| ). We thus obtain 
2/v^dif ~ 1-55, in agreement with the data of Figure [I], £dif ~ 1-61, i.e., 2/^£ di{ ps 1.57. 
We also obtain a ~ 1.5, albeit with a large uncertainty. 
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Figure 6: Logarithmic plot of the distribution p u of number of flips. Full symbols: data for 
v = 2k even. Open symbols: data for v = 2k + 1 odd. Statistical errors are much smaller than 
the symbol size. Full lines: common fit described in the text. 



Besides the above analysis of the regime of large numbers of spin flips, due to long sur- 
viving free spins, our data yield yet other interesting informations. First, the persistence 
probability, i.e., the probability for a spin to never flip, has the following value: 

p = 0.44739. 

Then, the probability that a spin flips an even or an odd number of times reads 
Pcvcn = E P^k = 0.72302, P odd = ]T p 2k+1 = 1 - P cvcn = 0.27698. 

fe>0 fc>0 

These figures can be related to the overlap between a random initial configuration and 
the corresponding final one. We have indeed 

Q = (<7 (0)<7 (oo)) = = Pcvc„ - P,dd = 0.44604. 

The mean number of spin flips reads 

(u) = Y J vVv = 2.06916. 

This quantity can be used to determine the fractions / CO nd of condensation moves and 
/dif of diffusive moves in a typical history, starting from a random initial configuration. 
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Indeed, on the one hand, a condensation move lowers the total energy by four units, 
while a diffusive moves leaves it unchanged (see Table [I]). On the other hand, any move 
involves exactly two spin flips. We have therefore, using 



Zcond = 1 - /dif = t7~t = 0.12477. 

2(u) 

To close up, let us compare the above results to the case of constrained Kawasaki 



dynamics (Wo = 0) || [L^l , where only condensation moves are allowed. In this situation, 



the above quantities can be simply related to the final energy E of (|4.1| ). Indeed every 



spin flips at most once, so that only po and p\ = 1 — p$ are non-zero, and f con d = 1- 
Starting again from a random configuration, we are left with 

2(l — p ) = l — Q = 2{v) = -E = 0.274087. 



6 Discussion 

First, we wish to emphasize the richness of the zero-temperature dynamics of the ferro- 
magnetic Ising chain. There are indeed four different natural kinds of dynamics, summa- 
rized in Table [| Only Glauber dynamics gives rise to a bona fide coarsening dynamics, 
obeying dynamical scaling with a typical domain size growing as L{t) ~ y/i. With the 
three other dynamics, the system is left in a metastable configuration after a relatively 
short blocking time. The constrained Glauber and Kawasaki dynamics have been ana- 
lyzed in our previous work ||. Each spin may flip at most once, before a global blocked 
state is reached. The blocking time scales with the number of spins as In TV. These 
models can be exactly mapped onto the random sequential addition problem of dimers 
and hollow trimers, respectively, hence allowing an analytical treatment. 



conserved 
magnetization 


diffusive 
moves 


dynamics 


behavior 


no 


yes 


Glauber |18| 












coarsening |19( 


no 


no 


constrained Glauber 
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met ast ability || 


yes 


yes 


Kawasaki |T 












metastability [this work] 


yes 


no 


constrained Kawasaki 


12 




met ast ability || 



Table 3: Four different zero-temperature dynamics of the ferromagnetic Ising chain, with 
appropriate references (brackets). 



The last case not considered so far from the viewpoint of metastability, Kawasaki 
dynamics, has been the subject of this work. Because of the diffusive moves, the dynamics 
is only partly irreversible. This novel feature makes the Kawasaki problem both richer 
and hopefully closer to more realistic situations. The number of flips of a given spin, 
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although finite with probability one, may be arbitrarily large. The blocking time grows 
as (In TV) 3 . On the other hand, analytical tools being no longer available, we needed to 
have recourse to numerical simulations. 

The present work demonstrates that there is no good a priori statistical description 
of the metastable states reached by Kawasaki dynamics. In other words, the Edwards 
hypothesis is invalidated, at least as an exact prescription. Systematic differences are 
indeed observed between numerical results and a priori predictions, especially in the 
pattern of spin correlations (Figure |j) and in the distribution of domain sizes (Figure |5|). 
The latter are neither statistically independent nor exponentially distributed. 

We also want to underline that the effective trapping description of the late stages 



of the zero-temperature Kawasaki dynamics, already emphasized in Ref. |[L1| , appears to 
yield quantitative predictions for several novel physical quantities, including the statistics 
of the blocking time (( |2.4| ), ( |2.5[ ), Figure 0), and the distribution of the number of spin 
flips ((§3), 0, Figure g). 
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